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ABSTRACT 



The increasing interest in long term (inter-annual) 
weather changes and their relation to processes in the ocean 
is beginning to illuminate the need for and the lack of long 
term records of physically significant variables which occur 
over the vast oceanic regions of the northern hemisphere. 

An attempt is made here to evaluate the accuracy of a hind 
cast time series of surface wind vector components and speeds 
for the period of January 1960 through December 1969. 

The quasi -geostrophi c model used to calculate these 
records is described as well as the twelve-hourly surface 
pressure data which were used as input. The central moments 
of the probability distributions of the computed records are 
compared to those of corresponding time series observed on 
Ocean Station Vessels. Linear correlation coefficients 
between observed and computed records were found to average 
0.81 for the components of the wind vector and 0.65 for the 
wind speed. Regression relations between computed and 
observed records also are presented. Spectral analysis of 
low pass filtered records showed coherence between observed 
and computed records increased with decreasing frequency. 

The accuracy of the computed records in low latitudes also is 
esti mated . 
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I INTRODUCTION 



A. GENERAL 

One of the most important variables encountered in the 
study of physical oceanography and air-sea interaction is 
the wind field at the sea surface. To a large degree, the 
surface wind velocity and its variability, both in time and 
space, determine the exchange of heat and momentum across 
the air-sea interface. Through such exchange the wind affects 
not only oceanic surface phenomena, such as waves and surface 
currents, but also features within the sea such as the 
vertical thermal structure and deep currents. 

The smallest scale of concern in this study is the 
synoptic scale. Such a scale covers a period of time from 
several hours to several days and a distance from several 
tens of kilometers to several hundred kilometers. Past 
attempts at estimating the surface wind field on this scale 
have been of two types. One could either average wind obser- 
vations collected over a given area for a specified period 
of time (Seckel, 1970) or one could estimate the surface winds 
from a similar collection of observations of other variables, 
such as surface atmospheric pressure (Roden, 1974). 

The accuracy of synoptic scale surface wind fields 
derived directly from wind observations suffers from two 
major drawbacks: a paucity of observations, and the high 
variability of the small scale wind field. The small number 
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of observations usually available in a typical mid-ocean 
region over a three-or-four-day period are very often grouped 
together in a small section of the region, often along major 
shipping lanes or along the track of a single ship. In 
addition, major portions of the open ocean lack any reports 
of wind conditions for weeks at a time. As a result, the 
few reports available must necessarily be used to describe 
large areas when estimating the synoptic scale wind field. 
Verploegh ( 1 967) found, however, that s i mul taneolis wind speed 
observations from ships less than sixty miles apart in 
"synopti cal ly homogeneous areas" differed significantly. His 
analysis indicated the correlation coefficient between such 
observations was less than 0.7. 

Because of such considerations, the second, indirect 
approach is generally used to estimate the synoptic scale 
surface wind field. This approach depends on the fact that 
the distribution of surface atmospheric pressure can be more 
accurately described than can the surface wind field from a 
limited number of observations, due to the lower variability 
in time and space of the surface atmospheric pressure field. 
An accurate estimate of the synoptic-scale surface wind field 
may be made by applying the geostrophic approximation to the 
surface atmospheric pressure field. 



B. 



OBJECTIVES 



Although the relation of the geostrophic to observed 
wind is a popular topic, this study is not intended to be a 
detailed examination of that relation. Rather, it is limited 
to an evaluation of one particular surface wind model used, 
in this study, to hindcast time series of surface winds over 
oceanic areas. The model is quasi -geostrophi c , i ncorporati ng 
the effects of friction. It is one part of operational 
oceanographic models used at Fleet Numerical Weather Central 
and at Fleet Weather Central, Rota (Spain). 

As a result of this study, the ten-year average and 
linear trend of observed wind conditions at several points 
in the North Atlantic and North Pacific Oceans were determined 
and compared to those calculated from the model. The primary 
objective of the study was however, to determine the accuracy 
of the computed time series and to determine if calibration 
factors could be applied directly to them in order to increase 
their accuracy. 



The model was developed by the author while at Fleet 
Numerical Weather Central; most of the work of this study was 
carried out while the author was at the Environmental Prediction 
Research Facility. 
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II. APPROACH 



A. THEORETICAL BACKGROUND 

1 . The Geostrophic Wind Equation 

Geostrophic wind has been considered to be an 
acceptable estimate of the true surface wind in many studies 
of oceanography and air-sea interaction (Roden 0£. c i t . , 
Namias, 1963). Geostrophy assumes a balance between the 
pressure gradient force and the Coriolis force and ignores 
acceleration, friction, and vertical motion. This balance is 
given in equation (1): 

<’) 

where = geostrophic wind vector 
P = surface pressure 

g g ~ 

V H = 3~x ^ + Ty (horizontal gradient operator) 

X = vector product operator 
f = 2 uj sin 0 (the Coriolis parameter) 
co = angular rotation of the earth 
0 = latitude 

p = air density 

A 

k = unit vertical vector 

Although this approximation is widely used, its accuracy is 
limited by the required assumptions indicated above. 
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The true surface wind vector is the resultant of 
the geostrophic wind vector and an ageostrophic wind vector. 

The direction and magnitude of the ageostrophic wind vector 
determines the accuracy of the geostrophic approximation. 
Unfortunately the ageostrophic vector has been shown to vary 
significantly over the Northern Hemisphere (Roll, 1965; 
Carstensen, 1967; Brummer, e_t aJL , 1 974). Such variability 
makes the geostrophic wind hard to relate to the surface wind 
in studies on synoptic or larger scales. 

The accuracy of the geostrophic approximation is 
also limited, in practical application, by the temporal 
variability of the surface pressure field. In practice, the 
geostrophic approximation is often applied to surface pressure 
distributions averaged over a month or more (Namias, op ci t . ) . 
Such a procedure results in calculating the vector resultant 
of the geostrophic wind for the averaging period. Obviously, 
non-conservative phenomena which depend on the surface wind 
cannot be treated using such a resultant wind field. In the 
study of such phenomena, the resultant vector is assumed to 
represent the mean or steady-state wind prevailing over the 
averagi ng interval . 

Whether or not the resultant wind is a good esti- 
mate of the mean wind depends on the variance of the wind 
vector over that interval. The mean wind speed, W^, and the 
resultant wind speed, W R are given below in terms of components 
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( 2 ) 



W 



M 



9 9 ^ 

t (ur + vf) 

i — - — 

t=l T 



W R 





(3) 



where U^, are orthogonal components of the wind vector at 

time t for 1 <_ t <_ T . (Averaging twelve hourly data over 

thirty days, we have 1 <_ t <_ 60). The difference between 
2 2 

and can be expressed in terms of the variance of each 

2 2 2 
of the components, a u and a y , and of the wind speed, 



,,2 ,.2 2.2 2 

W M ' W R °u + °v ' °w 



(4) 



where 



T 

Z 

t=l 



u t‘ 



/ 



u t \ 



t= 1 



) 



(5) 



2 2 

and similarly for a , and a . As will be shown in Section 
III -A , the difference expressed in equation (4) results in 
a drastic under-estimate of wind speed when the resultant 
wind is used in place of the mean wind speed in problems 
involving the square of the wind speed. 

2 . The Surface Wi nd-Geostrophi c Wind Relation 

Although this is not a study of the general relation 
of the geostrophic wind to the surface wind, a few words 
regarding this relation are appropriate for background. As 
indicated, the surface wind vector can be resolved into 
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geostrophic and ageostrophic components. The ageostrophic 
component can be assumed to consist of a frictional component 
and an acceleration component (Brummer, e_t a_l_. , op c i t . ) . 

The frictional component is generally accepted as 
the dominant ageostrophic component (Verploegh, op c i t . , 
Brummer, et a]_. , op c i t . ) . Haltiner and Martin (1 957 ) suggest 
that the frictional component is dependent on the stability 
of the surface layer of the atmosphere. Hasse (1974) has 
demonstrated that the dependence of the surface wind speed on 
the surface pressure gradient is an order of magnitude greater 
than its dependence on stability. 

The effects of stability have not been explicitly 
incorporated into the model; to some degree however, they 
have been implicitly incorporated through the latitude 
dependence of terms in the various equations of the model. 

The major effects of stability are assumed to occur in the 
region between 45°N and 25°N. This implicit incorporation of 
stability is shown in the ratio of the surface wind speed to 
the geostrophic wind speed in Figure 1. 

B. METHOD OF COMPUTING THE SURFACE WIND 

The calculations were performed on the Fleet Numerical 
Weather Central (FNWC) 63x63 polar stereographic grid of the 
northern hemisphere (FNWC, 1974). On this square grid the 
equator is an inscribed circle, i.e., tangent to the grid 
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boundary. The grid mesh length at 60°N is 381 km and increases 
with increasing latitude. The data, although calculated on a 

polar stereographic grid, will be discussed relative to a 

2 

Mercator grid. 

The components of the horizontal pressure gradient at a 
given grid point (i,j), were approximated to the second order 
by equation ( 6 ): 

fy(ij) = K n {P(i-2 , j ) - K 2 [P(i-l,j) 

( 6 ) 

- P ( i + 1 > j ) ] - P ( i + 2 , j ) } 

where = 1/12 

K 2 =8 

L = 381 km (mesh length) 

M F ( i , j ) = 1 0 ) — ~ = map factor of the grid at 

point (i ,j) 

0 = latitude of grid point (i,j), 

9 P 

and similarly for -ry(i,j). Using (1) with a constant air 

- 3 - 3 

density of 1.22 x 10 gm cm , the geostrophic wind vector 



2 

Both Mercator and polar stereographic projections are 
conformal projections of the earth and as a result, angles 
are preserved when transformed from one to the other (Taylor, 
1955) . 
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